Spin polarized states in nuclear matter 
with Skyrme effective interaction 



O 

o 



O 



A. A. Isayev 
Kharkov Institute of Physics and Technol 
Academicheskaya Str. 1, Kharkov, 61108, Ukraine 

J. Yang 

Dept. of Physics and Center for Space Science and Technology, 
Ewha Womans University, Seoul 120-750, Korea 
and Center for High Energy Physics, Kyungbook National University, Daegu 702-701, Korea 



The possibility of appearance of spin polarized states in symmetric and strongly asymmetric nu- 
clear matter is analyzed within the framework of a Fermi liquid theory with the Skyrme effective 
interaction. The zero temperature dependence of the neutron and proton spin polarization param- 
I eters as functions of density is found for SkM*, SGII (symmetric case) and SLy4, SLy5 (strongly 
■ asymmetric case) effective forces. By comparing free energy densities, it is shown that in symmetric 
nuclear matter ferromagnetic spin state (parallel orientation of neutron and proton spins) is more 
preferable than antiferromagnetic one (antiparallel orientation of spins). Strongly asymmetric nu- 
clear matter undergoes at some critical density a phase transition to the state with the oppositely 
directed spins of neutrons and protons while the state with the same direction of spins does not 
I , appear. In comparison with neutron matter, even small admixture of protons strongly decreases 
^ • the threshold density of spin instability. It is clarified that protons become totally polarized within 
' a very narrow density domain while the density profile of the neutron spin polarization parameter 
l^-^ is characterized by the appearance of long tails near the transition density. 

o 

^ PACS numbers: 21.65. +f; 75.25.+z; 71.10.Ay 

o ■ 

O ■ I. INTRODUCTION 



The spontaneous appearance of spin polarized states in nuclear matter is the topic of 
a great current interest due to relevance in astrophysics. In particular, the effects of spin 
^ ' correlations in the medium strongly influence the neutrino cross section and neutrino lu- 
(~| , minosity. Hence, depending on whether nuclear matter is spin polarized or not, drastically 
^ ' different scenarios of supernova explosion and cooling of neutron stars can be realized. An- 
. other aspect relates to pulsars, which are considered to be rapidly rotating neutron stars, 
surrounded by strong magnetic field. There is still no general consensus regarding the mech- 
?H [ anism to generate such strong magnetic field of a neutron star. One of the hypotheses is 
. . 1 that magnetic field can be produced by a spontaneous ordering of spins in the dense stellar 
core. 

The possibility of a phase transition of normal nuclear matter to the ferromagnetic state 
was studied by many authors. In the gas model of hard spheres, neutron matter becomes 
ferromagnetic at g ^ 0.41 fm~^ l]. It was found in Refs. 0, 01 that the inclusion of 
long-range attraction significantly increases the ferromagnetic transition density (e.g., up to 
g w 2.3 fm~'^ in the Brueckner theory with a simple central potential and hard core only for 
singlet spin states 3]). By determining magnetic susceptibility with Skyrme effective forces, 
it was shown in Ref. 4] that the ferromagnetic transition occurs at g 0.18-0.26 fm~^. The 
Fermi liquid criterion for the ferromagnetic instability in neutron matter with the Skyrme 
interaction is reached at p ~ 2-4po Q , where go = 0.16 fm~'^ is the nuclear matter saturation 
density. The general conditions on the parameters of neutron-neutron interaction, which 
result in a magnetically ordered state of neutron matter, were formulated in Ref. Spin 
correlations in dense neutron matter were studied within the relativistic Dirac-Hartree-Fock 
approach with the effective nucleon-meson Lagrangian in Ref. (3 , predicting the ferromag- 
netic transition at several times nuclear matter saturation density. The importance of the 
Fock exchange term in the relativistic mean-field approach for the occurrence of ferromag- 
netism in nuclear matter was established in Ref. 8]. The stability of strongly asymmetric 
nuclear matter with respect to spin fiuctuations was investigated in Ref. j9j, where it was 
shown that the system with localized protons can develop a spontaneous polarization, if the 
neutron-proton spin interaction exceeds some threshold value. This conclusion was con- 
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firmed also by calculations within the relativistic Dirac-Hartree-Fock approach to strongly 
asymmetric nuclear matter 10] . 

For the models with realistic nucleon-nucleon (NN) interaction, the ferromagnetic phase 
transition seems to be suppressed up to densities well above go 53-0. In particular, no 
evidence of ferromagnetic instability has been found in recent studies of neutron matter [T^ 
and asymmetric nuclear matter within the Brueckner-Hartree-Fock approximation with 
realistic Nijmegen II, Reid93 and Nijmegen NSC97e NN interactions. The same conclusion 
was obtained in Ref. , where magnetic susceptibility of neutron matter was calculated 
with the use of the Argonne vis two-body potential and Urbana IX three-body potential. 

Here we continue the study of spin polarizability of nuclear matter with the use of an 
effective NN interaction. As a framework of consideration, a Fermi liquid (FL) description 
of nuclear matter is chosen As a potential of NN interaction, we use the Skyrme 

effective interaction, utilized earlier in a number of contexts for nuclear matter calcula- 
tions ^3-j2i|. We explore the possibility of FM and AFM phase transitions in nuclear 
matter, when the spins of protons and neutrons are aligned in the same direction or in the 
opposite direction, respectively. In contrast to the approach, based on the calculation of 
magnetic susceptibility, we obtain the self-consistent equations for the FM and AFM spin 
order parameters and find their solutions at zero temperature. This allows us to determine 
not only the critical density of instability with respect to spin fluctuations, but also to 
establish the density dependence of the order parameters and to clarify the question of ther- 
modynamic stability of various phases. The main emphasis in our study will be laid on the 
region of zero isospin asymmetry (symmetric nuclear matter) and large isospin asymmetry 
(strongly asymmetric nuclear matter and neutron matter). 

Note that we consider the thermodynamic properties of spin polarized states in nuclear 
matter up to the high density region relevant for astrophysics. Nevertheless, we take into 
account the nucleon degrees of freedom only, although other degrees of freedom, such as 
pious, hyperons, kaons, or quarks could be important at such high densities. 



II. BASIC EQUATIONS 

The normal states of nuclear matter are described by the normal distribution function of 
nucleons /kik2 — Tr£)a+,aKj, where k = (p, cr, r), p is momentum, (7(t) is t he proi ection 
of spin (isospin) on the third axis, and g is the density matrix of the system [23l l24j. The 
energy of the system is specified as a functional of the distribution function f,E = E(f), 
and determines the single particle energy 

The self-consistent matrix equation for determining the distribution function / follows from 
the minimum condition of the thermodynamic potential and is 

/ = {exp{Yoe + Y^) + ly^ = {exp{YoO + 1}^' ■ (2) 

Here the quantities e and I4 are matrices in the space of k variables, with I4K1K2 = i^Ti<^KiK2 
(ti = n,j5), Yq = 1/T, Y4n = —^jf^/T and Yip = —^^/T being the Lagrange multipliers, 
/ij^ and /ip the chemical potentials of neutrons and protons, and T the temperature. Since 
it is assumed to consider a nuclear system with an excess of neutrons, the positive isospin 
projection is assigned to neutrons. Further we shall study the possibility of formation of 
various types of spin ordering (with parallel and antiparallel orientation of neutron and 
proton spins) in nuclear matter. 

The normal distribution function can be expanded in the Pauli matrices cr^ and in spin 
and isospin spaces 



/(p) = /oo(p)ctoTo + /3o(p)ct3To 

+ /o3(p)ctoT3 -I- /33(p)ct3T3. 



(3) 
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For the energy functional invariant with respect to rotations in spin and isospin spaces, the 
structure of the single particle energy is similar to the structure of the distribution function 
/: 

e(p) = eoo(p)croTo + e3o(p)CT3To (4) 

+ e03(p)f70T3 + e33(p)cr3T-3- 

Using Eqs. (j^J and (^J, one can express evidently the distribution functions /oo, /30, fo3, fas 
in terms of the quantities e: 

ho = ^{n{LJ+^+) + n{uj+-) + n{uj-^+) + ?i(tj_ _)}, 

/30 = ^{n{i^+.+) + -) - - 

fo3 = ^{n{i^+.+) - n{uj+-) + - 

/33 = ^{n{uj+^+) - _) - n{uj^^+) + _)}. (5) 

Here n{u!) — {exp(yot^) + 1}^^ and 

^+,+ = (.00 + £.30 + £.03 + £.33, 
^+,- = Coo + (,30 - (03 ~ (33, 
^-..+ = (00 — (30 + (03 — (33, 
~ (00 — (30 — (03 + (33, 



where 



(00 — £00 — Moo' ^30 — £30, 

(03 = £03 — M03' $33 = £33, 



,0 



t^OO — 2 ' ^^03 ^ 

As follows from the structure of the distribution functions /, the quantity w±,±, being the 
exponent in the Fermi distribution function n, plays the role of the quasiparticle spectrum. 
We consider the case when the spectrum is fourfold split due to the spin and isospin de- 
pendence of the single particle energy e(p) in Eq. Q). The branches <j-'±,+ correspond to 
neutrons with spin up and spin down, and the branches lu±,^ correspond to protons with 
spin up and spin down. 

The distribution functions / should satisfy the normalization conditions 

p 

^ 51 /o3(p) ^ Qn- Qp^ ag, (7) 
p 

^ .f3o(p) ^ - Qi^ , (8) 

p 

^ /33(p) = {gni + gpi) - {g,ii + jppt) = ^gn- (9) 
p 

Here a is the isospin asymmetry parameter, gn\, g-ni and Qp], Qpi are the neutron and proton 
number densities with spin up and spin down, respectively; p| = Pnt + (?pT ^'^'^ £*! ~ + 
are the nucleon densities with spin up and spin down. The quantities Ap|f and Ag|| may be 
regarded as FM and AFM spin order parameters. Indeed, in symmetric nuclear matter, if all 
nucleon spins are aligned in one direction (totally polarized FM spin state), then — g 

and A.gij^ = 0; if all neutron spins are aligned in one direction and all proton spins in the 
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opposite one (totally polarized AFM spin state), then Ap|j^ = q and A^?|-f = 0. In turn, 
from Eqs. one can find the neutron and proton number densities with spin up and 

spin down as functions of the total density isospin excess 5q = ag, and FM and AFM 
order parameters Agif and Ag||: 





1 , 


Sg + 




+ Apti), 




1 , 


5g — 


Aptt 


- Api^), 






5g + 


Apfi 


- Apii), 






5g — 


Apfi 


+ ApTi). 



In order to characterize spin ordering in the neutron and proton subsystems, it is conve- 
nient to introduce neutron and proton spin polarization parameters 

n„= np = M^^. (10) 

g-a gp 

The expressions for the spin order parameters Ag^^- and A^i^^ through the spin polarization 
parameters read 

AfillT = Sn^n + PpIIp, Ag-^i = gn^n - gp^p- 

To obtain the self-consistent equations, we specify the energy functional of the system in 
the form 

E{f)=Eo{f)+E„,t{f), (11) 

Eoif) - 45]eo(p)/oo(p), eo(p) = (12) 
p ° 

E.^tif) = 2^{eoo(p)/oo(p) + £3o(p)/3o(p) (13) 
p 

+ eo3(p)/o3(p) + e33(p)/33(p)}, 

eoo(p) ^^Yj ^o(k)/oo(q), k = 
q 

£30(p) = ^E^lW/30(q)' 

q 

^~03(p)-^E'^2(k)/03(q), 

q 

£33(p) = ^E^3W/33(q)- 

q 

Here toq is the bare mass of a nucleon, [/o(k), C/3(k) are the normal FL amplitudes, and 
£oOi £307 ^03) ^33 are the FL corrections to the free single particle spectrum. Further we do not 
take into account the effective tensor forces, which lead to coupling of the momentum and 
spin degrees of freedom p^L l2fil , and, correspondingly, to anisotropy in the momentum 
dependence of FL amplitudes with respect to the spin polarization axis. Using Eqs. 1^ and 
we get the self-consistent equations in the form 

Coo(p) = eo(p) + £oo(p) - Moo' 6o(p) = £3o(p), (14) 
^03 (p) = £03 (p) - ^J■o3^ 63 (p) = £33 (p)- 
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To obtain numerical results, we use the Skyrme effective interaction. In the case of Skyrme 
forces the normal FL amplitudes read 

C/o(k) = 6to + hQ^ + -|[3ii + i2(5 + Ax2)]k^, (15) 
h 



2 [ti(l - 2xi) - t2{l + 2x2)]\<? =a + b\f? 



[/2(k) = -2to(l + 2x0) - ^hQ^{l + 20:3) 



2[ti(l + 2a;i) -t2(l + 2a;2)]k2, 



C/3(k) = -2io - ^ise'^ - 4(^1 - i2)k2 = c + dk2, 

where ti,Xi,(3 are the phenomenological constants, characterizing a given parametrization 
of the Skyrme forces. Eqs. p5|l are derived in Appendix. In the numerical calculations 
we shall use SkM* j2^ and SGII potentials, designed for describing the properties of 
systems with small isospin asymmetry, and SLy4 and SLy5 potentials 30] , developed to fit 
the properties of systems with large isospin asymmetry. With account of the evident form 
of FL amplitudes and Eqs. (Q-® , one can obtain 

^00 = 7:, Moo, (16) 

2 

?03 = 7:, A^os, (17) 

6o = (a + 4)^ + 4<^'>^"' 

63 = (c + d^)^ + |(q^>33, (19) 

where the effective nucleon mass rrioo and effective isovector mass 77193 are defined by the 
formulae: 

■ -^[3ii+t2(5 + 4x2)], (20) 



2moo 27710 16 
^2 ^ 

Ty = ^[i2(l + 2x2)-ti(l + 2a;i)], 

277703 16 

and the renormalized chemical potentials /loo a-nd /103 should be determined from Eqs. © 
and 0. In Eqs. ifT^ and ifT^ . (q^)3o and (q^}33 are the second order moments of the 
corresponding distribution functions 

(q')30 = -^Eq'/3o(q), (21) 
q 

(q')33 = ^Eq'/33(q). (22) 



q 



Thus, with account of expressions (0 for the distribution functions /, we obtain the self- 
consistent equations Q^®, PT?! - and (|^ for the effective chemical potentials /^oo, M03, FM 
and AFM spin order parameters Agi||, Ap||, and second order moments (q^)3o, (q^)33- It is 
easy to see, that the self-consistent equations remain invariable under a global flip of spins, 
when neutrons (protons) with spin up and spin down are interchanged, and under a global 
flip of isospins, when neutrons and protons with the same spin projection arc interchanged. 

Let us consider, what differences will be in the case of neutron matter. Neutron matter is 
an infinite nuclear system, consisting of nucleons of one species, i.e., neutrons, and, hence. 
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the formalism of one-component Fermi liquid should be applied for the description of its 
properties. Formally neutron matter can be considered as the limiting case of asymmetric 
nuclear matter, corresponding to the isospin asymmetry a = 1. The individual state of 
a neutron is characterized by momentum p and spin projection a. The self-consistent 
equation has the form of Eq. (|2Jl, where all quantities are matrices in the space of k = (p, cr) 
variables. The normal distribution function and single particle energy can be expanded in 
the Pauli matrices in spin space 

/(p) = ./o(pko + /3(pV3, (23) 
£(p) = f^o(p)CTo +e3(p)CT3. 

The energy functional of neutron matter is characterized by two normal FL amplitudes 
C/(5'(k) and C/"(k). The normal FL amplitudes can be found in terms of the Skyrme force 
parameters ti,Xi,(3 (the details of the derivation procedure are given in Appendix): 

i7o"(k) = 2io(l-xo) + |/(l-X3) (24) 

+ ^[ti{l - Xl) + ^2(1 + X2)]k\ 

h 

U^{k) = -2toil-xo)-^-^g^il-x3) (25) 



2MI + X2) - ti(l - xi)]k2 = a„ + 6„k2 



With account of Eqs. (|24|l and (|25|l . the normalization conditions for the distribution func- 
tions can be written in the form 

p 

^5I■^3(p) = ^'T-a=A^?^■ (27) 

p 

Here g| and are the neutron number densities with spin up and spin down and 

f„ = ^{n{u;+)+n{u;^)}, c^±-eo±6, (28) 
/3 = i{n(^+) - n(a;_)}, (29) 
Co = -X Mn, (30) 

6 = («« + &„ — )^— + Y^(q )3- (31) 
The effective neutron mass to„ is defined by the formula 

^ = + f[Ml-.) + 3t.(l + ..)], (32) 

and the quantity (q^)3 in Eq. 131|) is the second order moment of the distribution function 
/a: 

(q')3 = |Eq'/3(q). (33) 
q 

Thus, with account of expressions (|28|l and (|29|l for the distribution functions /, we obtain 
the self-consistent equations (|?7jl . and (|^ for the effective chemical potential fin, spin 
order parameter Ag^-^, and second order moment (q^)3. 
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III. PHASE TRANSITIONS IN SYMMETRIC NUCLEAR MATTER 



Early research on spin polarizability of nuclear matter with the Skyrme effective inter- 
action were based on the calculation of magnetic susceptibility and finding its pole struc- 
ture determining the onset of instability with respect to spin fluctuations. Here we 
shall find directly solutions of the self-consistent equations for the FM and AFM spin order 
parameters as functions of density at zero temperature. In this section a special emphasis 
will be laid on the study of symmetric nuclear matter (a = 0), while in the next section on 
the investigation of strongly asymmetric nuclear matter {a ^ 1), including neutron matter 
as its limiting case. 

Let us consider the zero temperature behavior of spin polarization in symmetric nuclear 
matter {gp = Qn). The FM spin ordering corresponds to the case Ag^i ^ 0, (q^)3o 7^ 
0, A£i|| = 0, (q^)33 = 0, while the AFM spin ordering to the case Aq\i ^ 0, (q^)33 ^ 
0, Agf I — 0, (q^}3o = 0. In the totally ferromagnetically polarized state nontrivial solutions 
of the self-consistent equations have the form 

Aew = Q, (q')30 = ^Qkl. (34) 

Here kp = (Stt^^))^/"^ is Fermi momentum of symmetric nuclear matter in the case when 
degrees of freedom, corresponding to spin up of nucleons, are open while those related to 
spin down are inaccessible. For the totally antiferromagnetically polarized nuclear matter 




FIG. 1: FM and AFM spin polarization parameters as functions of density at zero temperature for 
(a) SkM* force and (b) SGII force. 
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we have 

^gn = 9: (q')33 = ^gkl. (35) 

The Fermi momentum kp is given by the same expression as in Eq. H34() since now de- 
grees of freedom, related to spin down of protons and spin up of neutrons, are inaccessible. 
The results of numerical determination of FM Ag^-i^/ g and AFM Agn/g spin polarization 
parameters are shown in Fig. for the SkM* and SGII effective forces. 

The FM spin order parameter arises at density g » 2go for the SkM* potential and at 
g « 2.75go for the SGII potential. The AFM order parameter originates at p « S.Spo for the 
SkM* force and a,t g ^ 5go for the SGII force. In both cases FM ordering appears earlier 
than AFM one. Nuclear matter becomes totally ferromagnetically polarized {Ag^-^/g — 1) 
at density g « 2.7go for the SkM* force and at p « 3.9go for the SGII force. Totally 
antiferromagnetically polarized state {Agn/g ~ 1) is formed at g « 4.5go for the SkM* 
potential and at g « 7.2po for the SGII potential. 

Note that in symmetric nuclear matter the neutron and proton spin polarization param- 
eters for the FM spin ordering are given by the formulas 

iiji — lip — 



and for the AFM spin ordering their expressions read 

_Agu 

llfl — lip — 



The second order moments (q^)3o, (q^)33 of the distribution functions /soi/aa also play 
the role of the order parameters. In Fig.[21it is shown behavior of these quantities normalized 
to their value in the totally polarized state. The ratios 5(q^)3o/3£ifc|, and 5{q^)33/3gkp are 
regarded as FM and AFM order parameters, respectively. The behavior of these quantities 
is similar to that of the spin polarization parameters in Fig. ^ with the same values of the 
threshold densities for the appearance and saturation of the order parameters. 

In the density domain, where FM and AFM solutions of self-consistent equations exist 
simultaneously, it is necessary to clarify, which solution is thermodynamically preferable. 
To this end, we should compare the free energies of both states. The results of the numerical 
calculation of the total energy per nucleon, measured from its value in the normal state, 
are shown in Fig. |3| One can see, that for all relevant densities FM spin ordering is more 
preferable than AFM one, and, moreover, the difference between corresponding free energies 
becomes larger with increasing density, so that there is no evidence, that AFM spin ordering 
could become preferable at larger densities. These results are in correspondence with the 
results of the Brueckner calculations with Nijmegen NSC97e potential in Ref. [l^j where it 
was shown that for symmetric nuclear matter the state with the oppositely directed spins 
of protons and neutrons is less favorable than the state with the same direction of nucleon 
spins for all relevant densities. However, in contrast to calculations in Ref. p^ . we find 
that for high density region there will be realized the FM spin ordering of nucleon spins as 
a ground state of symmetric nuclear matter instead of the nonpolarized state. 

IV. PHASE TRANSITIONS IN STRONGLY ASYMMETRIC NUCLEAR 

MATTER 

In this section we continue to study the properties of spin polarized nuclear matter, but 
now in the region of large isospin asymmetry. In contrast to symmetric nuclear matter, the 
analysis shows that in strongly asymmetric nuclear matter it is realized the antiparallel spin 
ordering of neutron and proton spins. 

If all neutron and proton spins are aligned in one direction, then for nontrivial solutions 
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FIG. 2: Same as in Fig. but for the second order moments of the distribution functions, 
of the self-consistent equations we have 



Qn — 
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(36) 



where kp = (Stt^p)^/'^ is the Fermi momentum of totaUy polarized symmetric nuclear matter. 
Therefore, for the partial number densities of nucleons with spin up and spin down one can 
get 



em: 



1 + a 



1 ~ a 



-9, PpT 



-Q, Qpl = Qui = 0. 



(37) 



If all neutron spins are aligned in one direction and all proton spins in the opposite one, 
then 



(38) 
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FIG. 3: Total energy per nucleon, measured from its value in the normal state, for FM and AFM 
spin ordering as a function of density at zero temperature for (a) SkM* force and (b) SGII force. 



and, hence, 

1 + a 1 — a 

PnT = ^ g, Qpl = ^ g, gpT = g"i 0- (39) 

Now we present the results of the numerical solution of the self-consistent equations with 
the effective SLy4 and SLy5 forces for strongly asymmetric nuclear {a = 0.95,0.9,0.8) and 
neutron (a = 1) matter. The neutron and proton spin polarization parameters n„ and Hp 
are shown in Fig. ^ as functions of density at zero temperature. Since in a polarized state 
the signs of spin polarizations are opposite, considering solutions correspond to the case, 
when spins of neutrons and protons are aligned in the opposite direction. Note that for SLy4 
and SLy5 forces, being relevant for the description of strongly asymmetric nuclear matter, 
there are no solutions corresponding to the same direction of neutron and proton spins. The 
reason is that the sign of the multiplier ts^—l + 2x3) in the density dependent term of the 
FL amplitude Ui, determining spin-spin correlations, is positive, and, hence, corresponding 
term increases with the increase of nuclear matter density, preventing instability with respect 
to spin fluctuations. Contrarily, the density dependent term —t^Q^/i m the FL amplitude 
C/3, describing spin-isospin correlations, is negative, leading to spin instability with the 
oppositely directed spins of neutrons and protons at higher densities. 

Another nontrivial feature relates to the density behavior of the spin polarization param- 
eters at large isospin asymmetry. As seen from Fig. 0] even small admixture of protons 
leads to the appearance of long tails in the density profiles of the neutron spin polarization 
parameter near the transition point to a spin ordered state. As a consequence, the spin 
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FIG. 4: Neutron and proton spin polarization parameters as functions of density at zero temperature 
for (a) SLy4 force and (b) SLy5 force. 

polarized state is formed much earlier in density than in pure neutron matter. For example, 
the critical density in neutron matter is g ~ 0.59 fm~^ for SLy4 potential and g w 0.58 fm~^ 
for SLy5 potential; in asymmetric nuclear matter with a = 0.95 the spin polarized state 
arises at g ~ 0.38 fm^^ for SLy4 and SLy5 potentials. Hence, even small quantity of protons 
strongly favors spin instability of highly asymmetric nuclear matter, leading to the appear- 
ance of states with the oppositely directed spins of neutrons and protons. As follows from 
Fig. ^ protons become totally spin polarized within a very narrow density domain (e.g., 
if a = 0.95, full polarization occurs at g ~ 0.41 fm^'^ for SLy4 force and at p « 0.40 fm^"^ 
for SLy5 force) while the threshold densities for the appearance and saturation of the neu- 
tron spin order parameter are substantially different (if a = 0.95, neutrons become totally 
polarized at g « 1.05 fm~'^ for SLy4 force and a,t g ^ 1.02 fm^'^ for SLy5 force). 
Note that the second order moments 

(q2)„ = (q2)„^ _ (q2)„^ (40) 
q 

q 

also characterize spin polarization of the neutron and proton subsystems. If the solutions 
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FIG. 5: Same as in Fig. 2] but for the second order moments (q^)p and (q^)n, normalized to their 
values in the totally polarized state. 

(q^}30 and (q^}33 of the self-consistent equations are known, then 

(q')n = ^((q'>30 + (q^)33), 

(q')p-^((q'>3o-(q^)33). 

The values of (q^)n and (q^)p for the totally polarized state are 

(q')no = ^gkUl + o.r\ (q')po = -Y^skUl af'^- 

In Fig.|Slwe plot the density dependence of the second order moments (q^}n and (q^}p, 
normalized to their values in the totally polarized state, for different asymmetries at zero 
temperature. These quantities behave similarly to the spin polarization parameters in Fig.0I 
i.e., there exist long tails in the density profiles of the neutron spin order parameter and the 
proton spin order parameter is saturated within a very narrow density interval. 

To check thermodynamic stability of the spin ordered state with the oppositely directed 
spins of neutrons and protons, it is necessary to compare the free energies of this state and 
the normal state. In Fig. |H| the difference between the total energies per nucleon of the 
spin ordered and normal states is shown as a function of density at zero temperature. One 
can see that nuclear matter undergoes a phase transition to the state with the oppositely 
directed spins of neutrons and protons at some critical density, depending on the isospin 
asymmetry. 
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FIG. 6; Total energy per nucleon, measured from its value in the normal state, for the state with 
the oppositely directed spins of neutrons and protons as a function of density at zero temperature 
for (a) SLy4 force and (b) SLy5 force. 

Note that our results of neutron matter calculations obtained with the Skyrme effective 
interaction predict a FM phase transition at some critical density, that is different from 
the results of calculations with realistic NN potentials jl^ J^]. In Ref. jl3|, employing 
Nijmegen II and Reid93 NN potentials, it has been found within the Brueckner-Hartree- 
Fock approach, that in the range of densities explored {g ^ 7go) totally polarized neutron 
matter is always more repulsive than nonpolarized matter, being an indication that a phase 
transition of neutron matter to a FM state is not expected. Besides, in Ref. jl^], in connec- 
tion with the problem of the neutrino diffusion in dense matter, it has been shown within 
the framework of the auxiliary field diffusion Monte Carlo method with the Argonne vig, 
two-body potential and Urbana IX three-body potential, that the magnetic susceptibility 
of neutron matter shows a strong reduction of about a factor 3 with respect to its Fermi 
gas value. However, calculations of the magnetic susceptibility with the Skyrme effective 
interaction '3| give indication of infinite discontinuity, and, hence, predict a FM transition 
at densities g ~ 0.18 — 0.26 fm~'^. 

For strongly asymmetric nuclear matter with Skyrme forces we find a phase transition to 
the state with antiparallel spins of neutrons and protons, that is different from the results 
of calculations with NSC97e NN potential in Ref. where the nonpolarized state was 
predicted for the whole range of densities up to 1.2 fm~'^. The reasons, explaining such 
disagreement in calculations with the effective and realistic NN potentials, are discussed in 
the next section. 
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V. DISCUSSION AND CONCLUSIONS 

Spin instability is a common feature, associated with a large class of Skyrme models, but 
is not realized in more microscopic calculations. The Skyrme interaction has been successful 
in describing nuclei and their excited states. In addition, various authors have exploited its 
applicability to describe bulk matter at densities of relevance to neutron stars 31]. The force 
parameters are determined empirically by calculating the ground state in the Hartree-Fock 
approximation and by fitting the observed ground state properties of nuclei and nuclear 
matter at the saturation density. In particular, the interaction parameters, describing spin- 
spin and spin-isospin correlations, are constrained from the data on isoscalar |32l l33j | and 
isovector (giant Gamow-Teller) |34ll35ll36l| spin-flip resonances. 

In a microscopic approach, one starts with the bare interaction and obtains an effective 
particle-hole interaction by solving iteratively the Bethe-Goldstone equation. In contrast 
to the Skyrme models, calculations with realistic NN potentials predict more repulsive total 
energy per particle for a polarized state comparing to the nonpolarized one for all relevant 
densities, and, hence, give no indication of a phase transition to a spin ordered state. It must 
be emphasized that the interaction in the spin- and isospin-dependent channels is a crucial 
ingredient in calculating spin properties of isospin symmetric and asymmetric nuclear matter 
and different behavior at high densities of the interaction amplitudes, describing spin-spin 
and spin-isospin correlations, lays behind this divergence in calculations with the effective 
and realistic potentials. 

In this study as a potential of NN interaction we chose SkM* and SGII (symmetric nuclear 
matter) as well as SLy4 and SLy5 (strongly asymmetric nuclear matter) Skyrme effective 
forces. The models SkM* and SGII 22j have been constrained by fitting the properties 
of nucleon systems with very small isospin asymmetries, while the models SLy4 and SLy5 
were further constrained to reproduce the results of microscopic neutron matter calculations 
(pressure versus density curve) (30||. Besides, in a recent publication it was shown that 
the density dependence of the nuclear symmetry energy, calculated up to densities q ^ 3po 
with SLy4 and SLy5 effective forces (together with some other sets of parameters among 
the total 87 Skyrme force parametrizations checked) gives the neutron star models in a 
broad agreement with the observables, such as the minimum rotation period, gravitational 
mass-radius relation, the binding energy, released in supernova collapse, etc. This is an 
important check for using these parametrizations in the high density region of strongly 
asymmetric nuclear matter. However, it is necessary to note, that the spin-dependent 
part of the Skyrme interaction at densities of relevance to neutron stars still remains to be 
constrained. Probably, these constraints will be obtained from the data on the time decay 
of magnetic field of isolated neutron stars 37]. In spite of this shortcoming, SLy4 and SLy5 
effective forces hold one of the most competing Skyrme parametrizations at present time 
for description of isospin asymmetric nuclear matter at high density (together with SkM* 
and SGII forces for description of symmetric nuclear matter) while a Fermi liquid approach 
with Skyrme effective forces provides a consistent and transparent framework for studying 
spin instabilities in a nucleon system. 

In summary, we have considered the possibility of phase transitions into spin ordered states 
of symmetric and strongly asymmetric nuclear matter within the Fermi liquid formalism, 
where NN interaction is described by the Skyrme effective forces (SkM*, SGII and SLy4, 
SLy5 potentials for the regions of vanishing and strong isospin asymmetry, respectively). 
In contrast to the previous considerations, where the possibility of formation of FM spin 
polarized states was studied on the base of calculation of magnetic susceptibility, we obtain 
the self-consistent equations for the FM and AFM spin order parameters and solve them in 
the case of zero temperature. It has been found that nuclear matter demonstrates different 
behavior at high densities with respect to spin fluctuations in isospin symmetric and strongly 
isospin asymmetric cases. In the model with SkM* and SGII effective forces symmetric 
nuclear matter undergoes a FM phase transition, when the spins of protons and neutrons are 
aligned along the same direction. In the model with SLy4 and SLy5 effective forces strongly 
asymmetric nuclear matter is subjected to a phase transition into the spin polarized state 
with the oppositely directed spins of neutrons and protons, while the state with the same 
direction of the neutron and proton spins docs not appear. In the last case, an important 
peculiarity of the corresponding phase transition is the existence of long tails in the density 
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profile of the neutron spin polarization parameter near the transition point. This means 
that even small admixtm'c of protons to neutron matter leads to a considerable shift of the 
critical density of spin instability in the direction of low densities. In the model with SLy4 
effective interaction this displacement is from the critical density g w S.Tpo for neutron 
matter io g k, 2AgQ for asymmetric nuclear matter at the isospin asymmetry a = 0.95, 
i.e. for 2.5% of protons only. As a result, the state with the oppositely directed spins of 
neutrons and protons appears, where protons become totally polarized in a very narrow 
density domain. 

As a consequence of this study, important questions appear, what is the value of the 
threshold asymmetry, at which the parallel spin ordering at small isospin asymmetry is 
changed to the antiparallel spin ordering at large isospin asymmetry, and do the obtained 
results survive for another type of an effective interaction, e.g., for Gogny effective force [sM 
I39I I or monopole effective interaction |4Q||? This research is in progress and will be reported 
elsewhere. 



Appendix 

The aim of this section is to establish relationships between the Fermi liquid amplitudes 
and amplitudes of NN interaction in the leading order on the interaction between nucleons. 
To this end we present the Hamiltonian of the system in the form 



H = Ho + V, (41) 
1 

2V 



where v is the amplitude of NN interaction. We shall assume that the amplitude v doesn't 
depend on the total momentum of colliding nucleons, but only from their relative momenta: 

?)(kiK2;k3'«4) = w(pi,P2;P3,P4)kiK2;k3K4, K = (cr, r), (42) 

w ^ \ P1-P2 P3-P4 

w(pi,P2;p3,P4)kiK2:k3k^ ^ v[p, <i)k,k^ik^k^; p = — 2 — ' — 2 — ■ 

To obtain the energy functional, corresponding to the Hamiltonian H41|) , we should average 
the Hamiltonian (|41|l over the state of nonideal gas of particles. In the leading approximation 
on the interaction, using the Wick rules and taking into account that Trfpaj^^aRj = /kik2 [q 
being the statistical operator), one gets 

E{f) ^ Eoif) + E.ntif) (43) 

^o(/) = ^ eKlK2/K2ftl7 Eintif) ^ — ^ ^'(kiK2; '«3K4)(/k3Ki/k4K2 ~ /k4Ki/k3K2)- 

Hence, according to Eq. {Q), expression for the single particle energy reads 

e«iK2 =e°4K2(/)+£«i«2(/), (44) 
£kiK2(/) = ^ ^{w(kiK3;k2K4) + -0(k3Ki; K4K2) - w(kiK3; K4K2) - w(k3Ki; K2K4)}/k4K3- 

Av3 Av4 

For spatially homogeneous distributions 

/kiK2 — /kiK2 (Pl)^Pl >P2 1 £kiK2 ~ (Pl )'5pi ,p2 (45) 

and Eqs. can be simphfied 

m)= E 44«2(p)42-i(p). (46) 

—\ '—2 ' P 



Eintif) — ^ E %iK2(P)'^«2«l(P) 
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where 



The general structure of the NN interaction ampUtude u, invariant with respect to rotations 
in spin and isospin spaces, has the form 

v{P,<i) = Wo(P,q) + wi(p,q)cricr2 +W2(p,q)TiT2 + t^sCP, q)(o-icr2)(TiX2). (48) 

Taking into account the general structure of the normal distribution function 

/(p) = /oo(p) + /feo(p)o-fe + /o((p)Ti + fki{p)(^kTi (49) 

and calculating traces in Eqs. (|46|l . 147(1 . one can get 

Eoif) - 4^eo(p)/oo(p), eo(p) = (50) 
p 

E^ntif) = 2^{£oo(p)/oo(p) + 4o(p)Ao(p) (51) 
p 

+ £oi{p)foi{p) + £ki{p)fki{p)}, 

eoo(p) = ^ H ^o(k)/oo(q), k = 
q 



£~fco(p) = ^E^i(k)Ao(q), 

q 

£oKp) = ^E^2(k)/o,(q), 
q 

£fci(p) ^ H ^3(k)/fci(q), 



2V 
q 



where 



{/o(k) = {4z;o(k, k) - «o(-k, k) - 3«i(-k, k) - 3i.2(-k, k) - 9v-,{-k, k)} + {k ^ -k}, 

C/i(k) = {4z;i(k,k) -wo(-k,k)+wi(-k,k) -3w2(-k,k)+3v3(-k,k)} + {k^ -k}, 

(52) 

C/2(k) = {4z;2(k, k) - z;o(-k, k) - 3z;i(-k, k) + v^{-k, k) + 3i>3(-k, k)} + {k ^ -k}, 
C/3(k) = {4z;3(k, k) - «o(-k, k) + vi{-k, k) + «2(-k, k) - i;3(-k, k)} + {k ^ -k}. 

The interaction energy functional l|13() represents a special case of the functional (|51|l . cor- 
responding to a coUinear spin ordering. The amplitude of NN interaction for the Skyrme 
effective forces has the form 

v{p, q) = toil + XoPa) + ^^3(1 + X3P„)qI^ + ^ii(l + 2:iP<x)(p' + q') (53) 

+ %{l + X2Pa)pC^, = i(l + (TiCTs), 

n ^ 

Here is the spin exchange operator. Extracting from Eq. H53|l the functions wq, and 
substituting their expressions into Eq. H52I) . we obtain Eqs. p5|) for the normal FL ampli- 
tudes Uq, ...,1/3, describing, respectively, density, spin, isospin and spin-isospin correlations 
in a nucleon Fermi liquid. 
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For neutron matter the basic equations l|43|l . (|44|l remain unchanged, but now the indi- 
vidual state of a neutron is characterized by momentum p and spin projection a, k = (p, cr). 
Besides, taking into account the general structure of the interaction amplitude v and the 
normal distribution function / 

z)(p,q) = Uo(p,q) +Wi(p,q)<Ti<T2, (54) 
/(p) = /o(p) + /fe(p)afc, (55) 

and calculating traces with the Pauli matrices (7^, we obtain 

,2 



Eoif) = 2^eo(p)/o(p), eo(p) - ^ , 
E^ntif) - ^{eo(p)/o(p) + £fc(p)/fe(p)}, 



where 



eoip) = ^Et^o"(k)/o(q), k= (56) 
q 

Skip) = ^Y.^l\k)Mci), 
q 

and the FL amplitudes are given by the formulas 

U^\-k) = {2i.o(k,k) -«o(-k,k) -3«i(-k,k)} + {k^ -k}, (57) 
U^ik) = {2t.i(k,k) ^ wo(-k,k) +?;i(-k,k)} + {k -> -k}. 

After substituting expressions for the functions vo,vi from Eq. (|53|l . we obtain Eqs. H24|l . 
(|25|l for the neutron matter FL amplitudes Uq, J7". 
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